A Brief History of Trans-Neptunian Space 



Eugene Chiang, Yoram Lithwick, and Ruth Murray-Clay 

University of California at Berkeley 

Marc Buie and Will Grundy 

Lowell Observatory 

Matthew Holman 

Harvard-Smithsonian Center for Astrophysics 

The Edgeworth-Kuiper belt encodes the dynamical history of the outer solar system. Kuiper belt objects 
(KBOs) bear witness to coagulation physics, the evolution of planetary orbits, and external perturbations 
from the solar neighborhood. We critically review the present-day belt's observed properties and the 
theories designed to explain them. Theories are organized according to a possible time-line of events. In 
chronological order, epochs described include (1) coagulation of KBOs in a dynamically cold disk, (2) 
formation of binary KBOs by fragmentary collisions and gravitational captures, (3) stirring of KBOs by 
Neptune-mass planets ("oligarchs"), (4) eviction of excess oligarchs, (5) continued stirring of KBOs by 
remaining planets whose orbits circularize by dynamical friction, (6) planetary migration and capture of 
Resonant KBOs, (7) creation of the inner Oort cloud by passing stars in an open stellar cluster, and (8) 
collisional comminution of the smallest KBOs. Recent work underscores how small, coUisional, primordial 
planetesimals having low velocity dispersion permit the rapid assembly of ~5 Neptune-mass oligarchs at 
distances of 15-25 AU. We explore the consequences of such a picture. We propose that Neptune-mass 
planets whose orbits cross into the Kuiper belt for up to ^-^20 Myr help generate the high-perihelion 
members of the hot Classical disk and Scattered belt. By contrast, raising perihelia by sweeping secular 
resonances during Neptune's migration might fill these reservoirs too inefficiently when account is made 
of how little primordial mass might reside in bodies having sizes on the order of 100 km. These and other 
frontier issues in trans-Neptunian space are discussed quantitatively. 



1. INTRODUCTION 

The discovery by Jewitt and Luu (1993) of what many 
now regard as the third Kuiper belt object opened a new 
frontier in planetary astrophysics: the direct study of trans- 
Neptunian space, that great expanse extending beyond the 
orbit of the last known giant planet in our solar system. 
This space is strewn with icy, rocky bodies having diam- 
eters ranging up to a few thousand km and occupying orbits 
of a formerly unimagined variety. 

Kuiper belt objects (KBOs) afford insight into processes 
that form and shape planetary systems. In contrast to main 
belt asteroids, the largest KBOs today have hfetimes against 
collisional disruption that well exceed the age of the uni- 
verse. Therefore their size spectrum may preserve a record, 
unweathered by erosive collisions, of the process by which 
planetesimals and planets coagulated. At the same time, 
KBOs can be considered test particles whose trajectories 
have been evolving for billions of years in a time-dependent 
gravitational potential. They provide intimate testimony of 
how the giant planets — and perhaps even planets that once 
resided within our system but have since been ejected — had 
their orbits sculpted. The richness of structure revealed by 
studies of our homegrown debris disk is unmatched by more 



distant, extra-solar analogues. 

Section 2 summarizes observed properties of the Kuiper 
belt. Some of the data and analyses concerning orbital ele- 
ments and spectral properties of KBOs are new and have not 
been published elsewhere. Section 3 is devoted to theoret- 
ical interpretation. Topics are treated in order of a possible 
chronology of events in the outer solar system. Parts of the 
story that remain missing or that are contradictory are iden- 
tified. Section 4 recapitulates a few of the bigger puzzles. 

Our review is packed with simple and hopefully illumi- 
nating order-of-magnitude calculations that readers are en- 
couraged to reproduce or challenge. Some of these con- 
firm claims made in the hterature that would otherwise find 
no support apart from numerical simulations. Many esti- 
mates are new, concerning all the ways in which Neptune- 
sized planets might have dynamically excited the Kuiper 
belt. While we outline many derivations, space limitations 
prevent us from spelling out all details. For additional guid- 
ance, see the pedagogical review of planet formation by 
Goldreich et al. (2004a, hereafter G04), from which our 
work draws liberally. 

2. THE KUIPER BELT OBSERVED TODAY 
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2.1. Dynamical Classes 

Outer solar system objects divide into dynamical classes 
based on how their trajectories evolve. Fig. 1 displays or- 
bital elements, time-averaged over 10 Myr in a numerical 
orbit integration that accounts for the masses of the four 
giant planets, of 529 objects. Dynamical classifications of 
these objects are secure according to criteria developed by 
the Deep Ecliptic Survey (DES) team (Elliot et al, 2005, 
hereafter £05). Fig. 2 provides a close-up view of a portion 
of the Kuiper Belt. We distinguish 4 classes: 

1. Resonant KBOs (122/529) exhibit one or more mean- 
motion commensurabilities with Neptune, as judged by 
steady libration of the appropriate resonance angle(s) {Chi- 
ang et al, 2003, hereafter C03). Resonances most heav- 
ily populated include the exterior 3:2 (Plutino), 2:1, and 
5:2; see Table 1. Of special interest is the first discov- 
ered Neptune Trojan (1:1). All Resonant KBOs (except the 
Trojan) are found to occupy e-type resonances; the abiUty 
of an e-type resonance to retain a KBO tends to increase 
with the KBO's eccentricity e (e.g., Murray and Dermott, 
1999). Unless otherwise stated, orbital elements are heUo- 
centric and referred to the invariable plane. Several (9/122) 
also inhabit inclination-type (i^) or mixed-type (ei^) reso- 
nances. None inhabits an eN-type resonance whose stability 
depends on the (small) eccentricity of Neptune. The latter 
observation is consistent with numerical experiments that 
suggest CN-type resonances are rendered unstable by adja- 
cent e-type resonances. 



Table 1 . — Observed Populations of Neptune Resonances 
(securely identified by the DES team as of 8 Oct 2005) 
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2. Centaurs (55/529) are non-Resonant objects whose 
perihelia penetrate inside the orbit of Neptune. Most Cen- 
taurs cross the Hill sphere of a planet within 10 Myr. Cen- 
taurs are Ukely descendants of the other 3 classes, recently 
dislodged from the Kuiper belt by planetary perturbations 
(Holman and Wisdom, 1993, hereafter HW9T); Tiscareno 
andMalhotra, 2003). They will not be discussed further. 

3. Classical KBOs (246/529) are non-Resonant, non- 
planet-crossing objects whose time-averaged (e) < 0.2 and 
whose time-averaged Tisserand parameters 



(T) = ((flN/a) + 2V(a/aN)(l - e^) cos Ai) (1) 

exceed 3. Here Ai is the mutual inclination between the or- 
bit planes of Neptune and the KBO, a is the semi-major axis 
of the KBO, and on is the semi-major axis of Neptune. In 
the circular, restricted, 3-body problem, test particles with 
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Fig. 2. — Same as Fig. 1, zoomed in. 

T > 3 and a > cannot cross the orbit of the planet 
(i.e., their perihelia q = a(l — e) remain greater than a^). 
Thus, Classical KBOs can be argued to have never under- 
gone close encounters with Neptune in its current nearly 
circular orbit and to be relatively pristine dynamically. In- 
deed, many Classical KBOs as identified by our scheme 
have low inclinations (i) < 5° ("cold Classicals"), though 
some do not ("hot Classicals"). Our defining threshold for 
(e) is arbitrary; like our threshold for (T), it is imposed 
to suggest — perhaps incorrectly — which KBOs might have 
formed and evolved in situ. 

Classical KBOs have spectral properties distinct from 
those of other dynamical classes: Their colors are more uni- 
formly red (Fig. 3; see the chapter by Cruikshank et al). 
According to the Kolmogorov-Smirnov test {Press et al, 
1992; Peixinho et al, 2004), the probabilities that Clas- 
sical KBOs have B — V colors and Boehnhardt-S' slopes 
{Boehnhardt et al, 2001) drawn from distributions identi- 
cal to those for Resonant KBOs are 10~^ and 10~^, re- 
spectively. When Classical KBOs are compared to Scat- 
tered KBOs (see below), the corresponding probabilities are 
10~^ and lO"**. An alternative interpretation is that low-i 
KBOs are redder than high-i KBOs {Trujillo and Brown, 
2002; Peixinho et al, 2004). This last claim is statistically 
significant when Classical, Scattered, and Resonant KBOs 
are combined and analyzed as one set (Fig. 4). However, no 
correlation between physical properties and i (or any other 
measure of excitation) has proven significant within any in- 
dividual class. 

Both the inner edge of the Classical belt at a w 37 AU, 
and the gap in the Classical belt at a w 40^2 AU and (i) < 
10° (see Fig. 2), reflect ongoing sculpting by the present- 
day planets. The iimer edge marks the distance out to which 
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Fig. 3. — Visual colors of KBOs and Centaurs calculated from 
published photometry, with the average uncertainty indicated by 
the upper left oval. The spectral slope S is calculated for wave- 
lengths in the range of Johnson V through /. Neutral (solar) colors 
are indicated by 0. Symbols for dynamical classes are the same as 
those for Fig. 1. Classical KBOs constitute a distinct red cluster, 
except for (35671), which also has a small semi-major axis of 38 
AU. Other classes are widely dispersed in color. 

the planets have eroded the Kuiper belt over the last few 
billion years (HW93; Duncan et al, 1995, hereafter D95). 
The gap is carved by the z^ig, vn, and secular resonances 
{HW92>; D95; Knezevic et al, 1991). At a secular resonance 
denoted by Vj, the orbital precession frequency of a test 
particle — apsidal if j < 10 and nodal if j > 10 — matches 
one of the precession eigenfrequencies of the planets (see 
chapter 6 of Murray and Dermott, 1999). For example, at 
low i, the f8 resonance drives e to Neptune-crossing values 
in ~10^yr. Particles having large i, however, can elude 
the 1^8 {Knezevic et al, 1991). Indeed, 18 KBOs of various 
classes and all having large (z) reside within the gap. 

By contrast, the outer edge of the Classical belt at a w 
47 AU is likely primordial. Numerous surveys (e.g., £05; 
Bernstein et al., 2004; and references therein) carried out 
after an edge was initially suspected (Jewitt et al., 1998) 
all failed to find a single object moving on a low-e or- 
bit outside 47 AU. The reality of the "Kuiper Cliff" is 
perhaps most convincingly demonstrated by Trujillo and 
Brown (2001), who simply plot the distribution of heliocen- 
tric discovery distances of (mostly Classical) KBOs after 
correcting for the bias against finding more distant, fainter 
objects. This distribution peaks at 44 AU and plummets to a 
value 10 times smaller at 50 AU. The statistical significance 
of the CUff hinges upon the fact that the bias changes less 
dramatically — only by a factor of 2.2-2.4 for reasonable pa- 
rameterizations of the size distribution — ^between 44 and 50 
AU. The possibility remains that predominantly small ob- 
jects having radii i? < 50 km reside beyond 47 AU, or that 
the Cliff marks the inner edge of an annular gap having ra- 
dial width > 30 AU {Trujillo and Brown, 2001). 

4. Scattered KBOs (106/529) comprise non-Classical, 



Fig. 4. — Spectral slope S vs. time-averaged inclination. The 
typical uncertainty in S is indicated by the dotted bar. Classical 
KBOs evince no trend of color with (t). The solid line is fitted to 
Classicals only; statistical tests using Spearman's rank-order coef- 
ficient and Kendall's tau (Press et al, 1992) show that no signif- 
icant correlation exists. When Classical, Resonant, and Scattered 
KBOs are combined, S and (i) correlate significantly (with a false 
alarm probability of 10^ '); the dashed line is a fit to all three 
classes. The two most neutral Classicals are (35671) and 1998 
WV24, having semi-major axes of 38 and 39 AU, respectively. 

non-Resonant objects whose perihelion distances q remain 
outside the orbit of Neptune. (The "Scattered-Near" and 
"Scattered-Extended" classes defined in EQ5 — see also 
Gladman et al. (2002) — are combined to simplify discus- 
sion. Also, while we do not formally introduce Oort cloud 
objects as a class, we make connections to that popula- 
tion throughout this review.) How were Scattered KBOs 
emplaced onto their highly elongated and inchned orbits? 
Appealing to perturbations exerted by the giant planets in 
their current orbital configuration is feasible only for some 
Scattered objects. A rule-of-thumb derived from numerical 
experiments for the extent of the planets' collective reach 
is 9 ^ 37 AU {D95; Gladman et al, 2002). Fig. 1 reveals 
that many Scattered objects possess q > 37 AU and are 
therefore problematic. Outstanding examples include 2000 
CR105 (q = 44 AU; Millis et al, 2002; Gladman et al, 
2002) and (90377) Sedna {q = 76 AU; Brown et al, 2004). 

These classifications are intended to sharpen analyses 
and initiate discussion. The danger lies in allowing them 
to unduly color our thinking about origins. For example, 
though Sedna is classified above as a Scattered KBO, the 
history of its orbit may be distinct from those of other Scat- 
tered KBOs. We make this distinction explicit below. 

2.2. Sky Density and Mass 

We provide estimates for the masses of the Kuiper belt 
(comprising objects having q ^ 60 AU and a > 30 AU; 
§2.2.1); the inner Oort Cloud (composed of Sedna-hke ob- 
jects; §2.2.2); and Neptune Trojans (a w 30 AU; §2.2.3). 
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2.2.1. Main Kuiper Belt 

Bernstein et al. (2004, hereafter BQA) compile data from 
published surveys in addition to their own unprecedent- 
edly deep Hubble Space Telescope (HST) survey to com- 
pute the cumulative sky density of KBOs versus appar- 
ent red magnitude tur ("luminosity function"), shown in 
Fig. 5. Sky densities are evaluated near the ecliptic plane. 
Objects are divided into two groups: "CKBOs" (similar 
to our Classical population) having heliocentric distances 
38 AU < d < 55 AU and ecliptic inclinations i < 5°, and 
"Excited" KBOs (similar to our combined Resonant and 
Scattered classes) having 25 AU < < 60 AU and i > 5°. 
Given these definitions, their analysis excludes objects with 
ultra-high perihelia such as Sedna. With 96% confidence, 
B04 determine that CKBOs and Excited KBOs have differ- 
ent luminosity functions. Moreover, neither function con- 
forms to a single power law from niR = 18 to 29; in- 
stead, each is well-fitted by a double power law that flat- 
tens towards fainter magnitudes. The flattening occurs near 
mR « 24 for both groups. To the extent that all objects 
in a group have the same albedo and are currently located 
at the same heliocentric distance, the luminosity function is 
equivalent to the size distribution. We define q as the slope 
of the differential size distribution, where dN oc R~^dR 
equals the number of objects having radii between R and 
R + dR. As judged from Fig. 5, for CKBOs, q flattens 
from 5.5-7.7 (95% confidence interval) to 1.8-2.8 as R de- 
creases. For Excited KBOs, q flattens from 4.0-4.6 to 1.0- 
3.1. Most large objects are Excited (see also Fig. 1). 

By integrating the luminosity function over all magni- 
tudes, B04 estimate the total mass in CKBOs to be 




where all CKBOs are assumed to have the same albedo p, 
heliocentric distance d, and internal density p. The solid 
angle subtended by the belt of CKBOs is A. Given uncer- 
tainties in the scaUng variables — principally p and p (see 
the chapter by Cmikshank et al. for recent estimates) — this 
mass is good to within a factor of several. The mass is con- 
centrated in objects having radii i? 50 km, near the break 
in the luminosity function. 

The mass in Excited KBOs cannot be as reliably cal- 
culated. This is because the sample is heterogeneous — 
comprising both Scattered and Resonant KBOs having a 
wide dispersion in d — and because corrections need to be 
made for the observational bias against finding objects near 
the aphelia of their eccentric orbits. The latter bias can be 
crudely quantified as (Q/q)'^^^: the ratio of time an ob- 
ject spends near its apheUon distance Q (where it is unde- 
tectable) versus its perihelion distance q (where it is usually 
discovered). An order-of-magnitude estimate that accounts 
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Fig. 5. — Cumulative sky density vs. apparent red magnitude for 
CKBOs and Excited KBOs, from B04. If N{< niR) oc lO"™", 
then the size distribution index g = 5q + 1 (see §2.2. 1). Envelopes 
enclose 95% confidence intervals. The top abscissa is modified 
from B04; here it assumes a visual albedo of 10%. Figure provided 
by Gary Bernstein. 

for how much larger A, d^, and (Q/q)^^"^ are for Excited 
KBOs than for CKBOs suggests that the former population 
might weigh ~ 10 times as much as the latter, or ~O.O5M0. 
This estimate assumes that q for Excited KBOs is such that 
most of the mass is concentrated near niji w 24, as is the 
case for CKBOs. If q for the largest Excited KBOs is as 
smaU as 4, then our mass estimate increases by a logarithm 
to ~O.15M0. 

2.2.2. Inner Oort Cloud (Sedna-Like Objects) 

What about objects with unusually high perihelia such as 
Sedna, whose mass is Ms ~ 6 x 10-''(i?/750km)3Me? 
The Caltech Palomar Survey searched / ^ 1 /5 of the celes- 
tial sphere to discover one such object (Brown etal, 2004). 
By assuming Sedna-like objects are distributed isotropically 
(not a forgone conclusion; see §3.7), we derive an upper 
limit to their total mass of Ms{Q / qf^ f~'^ ~ O.IM®. If 
all objects on Sedna-hke orbits obey a size spectrum resem- 
bling that of Excited KBOs (B04), then we revise the upper 
limit to ~O.3M0. The latter value is 20 times smaller than 
the estimate by Brown et al. (2004); the difference arises 
from our use of a more reaUstic size distribution. 

2.2.3. Neptune Trojans 

The first Neptune Trojan, 2001 QR322 (hereafter "QR"), 
was discovered by the DES (C03). The distribution of DES 
search fields on the sky, coupled with theoretical maps of 
the sky density of Neptune Trojans (Nesvorny and Dones, 
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2002), indicate that N ~ 10-30 objects resembling QR li- 
brate on tadpole orbits about Neptune's forward Lagrange 
(L4) point (O03). Presumably a similar population exists at 
L5. An assumed albedo of \2-A% yields a radius for QR of 
i? ~ 65-1 15 km. Spreading the inferred population of QR- 
like objects over the area swept out by tadpole orbits gives a 
surface mass density in a single Neptune Trojan cloud that 
approaches that of the main Kuiper beh to within factors of 
a few {Chiang andLithwick, 2005). Large Neptune Trojans 
are at least comparable in number to large Jupiter Trojans 
and may outnumber them by a factor of -^10 (CDS). 

2.3. Binarity 

Veillet et al. (2002) optically resolved the first binary 
(1998 WW31) among KBOs having sizes R ~ 100 km. 
Over twenty binaries having components in this size range 
are now resolved. Components typically have comparable 
brightnesses and are separated (in projection) by 300-10^ 
km (e.g., Stephens and Noll, 2006). These properties re- 
flect observational biases against resolving binaries that are 
separated by < 0.1" and that contain faint secondaries. 

Despite these selection biases, Stephens and Noll (2006) 
resolved as many as 9 out of 81 KBOs (~10%) with HST. 
They further report that the incidence of binarity appears 
4 times higher in the Classical disk than in other dynami- 
cal populations. It is surprising that so many binaries exist 
with components widely separated and comparably sized. 
A typical binary in the asteroid belt, by contrast, comprises 
dissimilar masses separated by distances only slightly larger 
than the primary's radius. Another peculiarity of KBO bi- 
naries is that components orbit one other with eccentricities 
of order unity. In addition, binary orbits are inchned rela- 
tive to their heliocentric orbits at seemingly random angles. 
See Noll (2003) for more quantitative details. 

Although close binaries cannot be resolved, their com- 
ponents can eclipse each other. Sheppardand Jewitt (2004) 
highlight a system whose light curve varies with large am- 
plitude and httle variation in color — suggesting that it is a 
near-contact binary. They infer that at least 10% of KBOs 
are members of similarly close binaries. 

Among the four largest KBOs having R iv 1000 km — 
2003 UB313, Pluto, 2005 FY9, and 2003 ELei— three are 
known to possess satellites. Secondaries for 2003 ELgi and 
2003 UB313 are 5% and 2% as bright as their primaries, and 
separated by 49500 and 36000 km, respectively (Brown et 
ah, 2005b). In addition to harboring Charon (Christy and 
Harrington, 1978), Pluto possesses two, more distant com- 
panions having i? ~ 50 km (Weaver et al. , 2006). The three 
satellites' orbits are nearly co-planar; their semi-major axes 
are about 19600, 48700, and 64800 km; and their eccentric- 
ities are less than 1% (Buie et al, 2006). 

3. THEORETICAL TIMELINE 

We now recount a possible history of trans-Neptunian 
space. Throughout our narration, it is helpful to remember 



that the timescale for an object of size R and mass M to 
colUde into its own mass in smaller objects is 
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where M is the rate at which mass from the surrounding 
disk impacts the object, a is the disk's surface mass density 
(mass per unit face-on area) in smaller objects, Urei is the 
relative collision velocity, h ^ Wrd/^^ is the effective verti- 
cal scale height occupied by colhders, and Q is the orbital 
angular frequency. Relative velocities t;rei depend on how 
e and i are distributed. Equation (3) requires that e's and 
i's be comparably distributed and large enough that gravi- 
tational focussing is ignorable. While these conditions are 
largely met by currently observed KBOs, they were not dur- 
ing the primordial era. Our expressions below represent ap- 
propriate modifications of (3). 

3.1. Coagulation 

The mass inferred for the present-day Kuiper belt, 
--0.05-0.3 Af® (§2.2), is wefl below that thought to have 
been present while KBOs coagulated. Kenyon and Luu 
(1998, 1999) and Kenyon (2002), in a series of particle- 
in-a-box accretion simulations, find that ^3-30 of 
primordial solids, spread over an annulus extending from 
32 to 38 AU, are required to (a) coagulate at least 1 ob- 
ject as large as Pluto and (b) coagulate ~10^ objects having 
i? > 50 km. The required initial surface density, a ~ 0.06- 
0.6gcm~^, is of order that of the condensible portion 
of the minimum-mass solar nebula (MMSN) at 35 AU: 
o-MMSN ~ 0.2gcm-2. 

3.1.1. The Missing-Mass Problem 

That primordial and present-day masses differ by 2 or- 
ders of magnitude is referred to as the "missing-mass" prob- 
lem. The same accretion simulations point to a possible 
resolution: Only ~ 1-2% of the primordial mass accretes to 
sizes exceeding ~100 km. The remainder stalls at comet- 
like sizes of ^0.1-10 km. Stunting of accretion is attributed 
to the formation of several Pluto-sized objects whose grav- 
ity amplifies velocity dispersions so much that collisions 
between planetesimals are erosive rather than accretionary. 
Thus, accretion in the Kuiper belt may be self-hmiting 
(Kenyon and Luu, 1999). The bulk of the primordial mass, 
stalled at cometary sizes, is assumed by these authors to 
erode away by destructive collisions over Gyr timescales. 

We can verify analytically some of Kenyon and Luu's 
results by exercising the "two-groups method" (G04), 
whereby the spectrum of planetesimal masses is approx- 
imated as bimodal. "Big" bodies each of size R, mass M, 
Hill radius i?H> and surface escape velocity fcsc comprise 
a disk of surface density S. They are held primarily re- 
sponsible for stirring and accreting "small" bodies of size 
s, surface density a, and random velocity dispersion u. By 
random velocity we mean the non-circular or non-planar 
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component of the orbital velocity. As such, u is propor- 
tional to the root-mean-squared dispersion in e and i. 
To grow a big body takes time 



R 

R 



pR 



(4) 



where the term in parentheses is the usual gravitational fo- 
cussing factor (assumed < 1). Gravitational stirring of 
small bodies by big ones balances damping of relative ve- 
locities by inelastic colhsions amongst small bodies. This 
balance sets the equilibrium velocity dispersion u: 
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Combining (4) and (5) implies 
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at a distance of 35 AU. All bodies reside in a remark- 
ably thin, dynamically cold disk: Eccentricities and inch- 
nations are at most of order u/{ila) ^ 0.001. Our nominal 
choices for a, E, and s are informed by Kenyan and Lmm's 
proposed solution to the missing-mass problem. Had we 
chosen values resembling those of the Kuiper belt today — 
X; ~ (7 ~ O.OIctmmsn — coagulation times would exceed 
the age of the solar system. 

The above framework for understanding the missing- 
mass problem, while promising, requires development. 
First, account needs to be made for how the formation 
of Neptune — and possibly other planet-sized bodies — 
influences the coagulation of KBOs. None of the sim- 
ulations cited above succeeds in producing Neptune-mass 
objects. Yet minimum-mass disks may be capable, in theory 
if not yet in simulation, of producing several planets having 
masses approaching that of Neptune at distances of 15-25 
AU on timescales much shorter than the age of the solar 
system (G04; see §§3.4-3.5). The inability of simulations 
to produce ice giants may arise from their neglect of small- 
s, low-u particles that can be efficiently accreted (G04). 
Sizes s as small as centimeters seem possible. How would 
their inclusion, and the consequent formation of Neptune- 
mass planets near the Kuiper belt, change our understand- 
ing of the missing-mass problem? Second, how does the 
outer solar system shed ~99% of its primordial sohds? The 
missing-mass problem translates to a "clean-up" problem, 
the solution to which will involve some as yet unknown 
combination of colhsional comminution, diffusive transport 
by interparticle collisions, gravitational ejection by planets, 
and removal by gas and/or radiation drag. 

3.1.2. The Outer Edge of the Primordial Planetesimal Disk 



How far from the Sun did planetesimals coagulate? The 
outer edge of the Classical disk at 47 AU (§2) suggests that 
planetesimals failed to form outside this distance. Extra- 
solar disks are also observed to have well-defined bound- 
aries. The debris disks encircling /3 Pictoris and AU Mi- 
croscopii exhibit distinct changes ("breaks") in the slopes 
of their surface brightness profiles at stellocentric distances 
of 100 AU and 43 AU, respectively (e.g., Kalas and Jewitt, 
1995; KrLst et ai, 2005). This behavior can be explained 
by having dust-producing parent bodies reside only at dis- 
tances interior to break radii {Strubbe and Chiang, 2006). 

We cannot predict with confidence how planetesimal 
disks truncate. Our understanding of how micron-sized 
dust assembles into the "small," super-meter-sized bodies 
that coagulation calculations presume as input is too poor. 
Recent work discusses how sohd particles might drain to- 
wards their central stars by gas drag, and how the accu- 
mulation of such solids at small stellocentric distances trig- 
gers self-gravitational collapse and the formation of larger 
bodies (Youdin and Shu, 2002; Youdin and Chiang, 2004). 
These ideas promise to explain why planetesimal disks have 
sharp outer edges, but are subject to uncertainties regarding 
the viability of gravitational instability in a turbulent gas. 
To sample progress on planetesimal formation, see Garaud 
and Lin (2004), Youdin and Goodman (2005), and Gomez 
and Ostriker (2005). In what follows, we assume that ob- 
jects having R ~ 100 km coagulated only inside 47 AU. 

3.2. Formation of Binaries 

To have formed from a fragmentary collision, binary 
components observed today cannot have too much angular 
momentum. Consider two big bodies undergoing a gravita- 
tionally focussed collision. Each body has radius R, mass 
Af, and surface escape velocity v^sc- Prior to the collision, 
their angular momentum is at most imax ~ Mv^scR- Af- 
ter the collision, the resultant binary must have angular mo- 
mentum L < Ln^ii^. Unless significant mass is lost from the 
colhsion, components can be comparably sized only if their 
separation is comparable to their radii. Pluto and Charon 
meet this constraint. Their mass ratio is ^1/10, their sep- 
aration is ~20i?piuto> and hence their angular momentum 
obeys L/Lj^ax ~ \/20/10 < 1. Canup (2005) explains 
how Charon might have formed by a collision. The re- 
maining satelhtes of Pluto {Stern et al., 2006), the satel- 
lites of Pluto-sized KBOs 2003 ELei and 2003 UB313, and 
the candidate near-contact binaries discovered by Sheppard 
and Jewitt (2004) might also have formed by collisions. 

By contrast, binary components having wide separations 
and comparable masses have too much angular momentum 
to have formed by gravitationally focussed collisions. And 
if collisions were unfocussed, collision times would exceed 
the age of the solar system — assuming, as we do throughout 
this review, that the surface density of big (i? ~ 100 km) 
bodies was the same then as now (§1; §3.9). 

Big bodies can instead become bound ("fuse") by purely 
gravitational means while they are still dynamically cold. 
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Indeed, such binaries testify powerfully to the cold state of 
the primordial disk. To derive our expressions below, recall 
that binary components separated by the Hill radius ^Ru 
orbit each other with the same period that the binary's cen- 
ter of mass orbits the Sun, ~f2~^. Furthermore, we assume 
that the velocity dispersion v of big bodies is less than their 
HiU velocity L'H = ili?H- Then big bodies undergo runaway 
cooUng by dynamical friction with small bodies and settle 
into an effectively two-dimensional disk (G04). Reaction 
rates between big bodies must be calculated in a 2-D geom- 
etry. Because u > v-a, reaction rates involving small bodies 
take their usual forms appropriate for three dimensions. 

Goldreich et al. (2002, hereafter G02) describe two col- 
lisionless formation scenarios, dubbed and LP's. Both 
begin when one big body (L) enters a second big body's 
(L) Hill sphere. Per big body, the entry rate is 



Rh 
R 



(8) 



'pR \ R J pR 

where a = R/Rh « 1.5 x 10-*(35 AU/o). If no other 
body participates in the interaction, the two big bodies 
would pass through their Hill spheres in a time ~rj~^ (as- 
suming they do not collide). The two bodies fuse if they 
transfer enough energy to other participants during the en- 
counter. In L^, transfer is to a third big body: L + L + L ^ 
-|- L. To just bind the original pair, the third body must 
come within Rn of the pair. The probability for this to 
happen in time fi"^ is P^a ~ Nu^~^. If the third body 
succeeds in approaching this close, the probabiUty that two 
bodies fuse is order unity. Therefore the timescale for a 
given big body to fuse to another by is 



1,1,3 



pR 



2 4 

a 



n 



2Myr, (9) 



where the numerical estimate assumes R = 100 km, S = 
O.OIctmmsn. and a = 35 AU. 

In L^s, energy transfer is to small bodies by dynami- 
cal friction: L + L + s°° L"^ + s'^. In time the 
pair of big bodies undergoing an encounter lose a fraction 
{ai}/ pR){vesc/uy^~^ of their energy, under the assump- 
tion Vesc > u > vn (G04). This fraction is of order the 
probability Pf,2g that they fuse, whence 



flPNbin I K^iiH ' to conclude that the steady-state fraction of 
KBOs that are binaries with separation Ru is 



/bin (a; ^ Rk) 



bin 



all 



pR 



0.4%. (11) 



As X decreases, shrinkage slows. Therefore /bin increases 
with decreasing x. Scaling relations can be derived by argu- 
ments similar to those above. If L'^s dominates, /bin oc 
for X > Ru{vh/u)'^ and /bin oc x~^ for x < Ru{vh/u)^ 
(G02). If dominates, /bin oc x"^/^. 

Alternative formation scenarios require, in addition 
to gravitational scatterings, physical collisions. Weiden- 
schilling (2002) suggests a variant of in which the third 
big body collides with one member of the scattering pair. 
Since physical collisions have smaller cross-sections than 
gravitational interactions, this mechanism requires -^10^ 
more big (R ^ 100 km) bodies than are currently observed 
to produce the same rate that is cited above for (Weiden- 
schilling, 2002). Funato et al. (2004) propose that observed 
binaries form by the exchange reaction Ls + L ^ L"^ + s: 
A small body of mass to, originally orbiting a big body of 
mass M, is ejected by a second big body. In the major- 
ity of ejections, the small body's energy increases by its 
orbital binding energy leaving the big bodies 

bound to each other with separation x {M/m)R. The 
rate-limiting step is the formation of the pre-existing (Ls) 
binary, which requires (as in the asteroid belt) two big bod- 
ies to collide and fragment. Hence 



pR 



tfuse.exchange ~ ^"^^^ 0.6 Myr. 



(12) 



Estimating /bin as a function of x under the exchange 
hypothesis requires knowing the distribution of fragment 
masses to. Whether L^s, L^, or exchange reactions domi- 
nate depends on the uncertain parameters E, a, and s. 

As depicted above, newly formed binaries should be 
nearly co-planar with the 2-D disk of big bodies, i.e., bi- 
nary orbit normals should be nearly parallel. Observations 
contradict this picture (§2.3). How dynamical stirring of the 
Kuiper belt subsequent to binary formation affects binary 
inclinations and eccentricities has not been investigated. 

3.3. Early Stirring by Growing Planetary Oligarchs 



1 



.f^y^^^^rMyr, (10) 



where we have used (5) and set s = 100 m. 

Having formed with semi-major axis x ^ Rh ^ 
7000i?, a binary's orbit shrinks by further energy trans- 
fer. If is the more efficient formation process, pass- 
ing big bodies predominantly harden the binary; if L'^s 
is more efficient, dynamical friction dominates hardening. 
The probability P per orbit that x shrinks from '^Ru to 
^Ryi/2 is of order either Pj^s or PL2g. We equate the for- 
mation rate of binaries, N^a/t fuse, with the shrinkage rate. 



Coagulation of KBOs and fusing of binaries cannot pro- 
ceed today, in part because velocity dispersions are now 
so large that gravitational focussing is defeated on a wide 
range of length scales. What stirred the Kuiper belt? There 
is no shortage of proposed answers. Much of the remaining 
review (§§3.3-3.7) explores the multitude of non-exclusive 
possibilities. We focus on stirring "large" KBOs Uke those 
currently observed, having R ^ 100 km. Our setting re- 
mains the primordial disk, of whose mass large KBOs con- 
stitute only a small fraction (1-2%; §3.1.1). 

Neptune and Uranus are thought to accrete as oligarchs, 
each dominating their own annulus of full width ~5 Hill 
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radii (G04; Ida and Makino, 1993; Greenberg et ah, 1991; 
the coefficient of 5 presumes that oligarchs feed in a shear- 
dominated disk in which planetesimals have random ve- 
locities u that are less than the oligarch's Hill velocity 
vh = rii?H- If u > ''H, oligarchs' feeding annuli are wider 
by ^u/vu- In practice, u/vh does not greatly exceed unity 
since it scales weakly with input parameters.) Each oligarch 
grows until its mass equals the isolation mass, 

Mp ~ 2na x 5Rh,p x a , (13) 

where -Rh.p is the oligarch's Hill radius. For o = 25 AU 
and Mp equal to Neptune's mass Mn = 17M(^, equa- 
tion (13) implies a ~ 0.9gcm~^ ~ Sctmmsn- About 5 
Neptune-mass oligarchs can form in nested annuli between 
15 and 25 AU. Inspired by G04 who point out the ease with 
which ice giants coagulate when the bulk of the disk mass 
comprises very small particles (§3.1.1), we assume that all 5 
do form in a disk that is a few times more massive than the 
MMSN and explore the consequences of such an initially 
packed system. 

While oligarchs grow, they stir large KBOs in their 
immediate vicinity. A KBO that comes within distance 
b of mass Mp has its random velocity excited to vk ^ 
(GMp/&)^/^. Take the surface density of perturbers to 
be Ep. Over time t, a KBO comes within distance b ~ 
[Mp/(Spf2i)]V2 of aperturber. Therefore 

VK ~ G^/^(MpEpm)i/'^ . (14) 

Since Neptune and Uranus contain more hydrogen than 
can be explained by accretion of icy solids alone, they 
must complete their growth within tacc,p ^ 

1-10 Myr, 

before all hydrogen gas in the MMSN photo-evaporates 
(e.g., Matsuyama et al, 2003, and references therein). For 
t = taccp = 10 Myr, Mp = Mn, Sp = 0.9 gcm-^, and 
fi = 27r/(100yr), equation (14) implies vk ~ lkms~^ 
or ex ~ 0.2. It is safe to neglect damping of vk for large 
KBOs, which occurs by inelastic colUsions over a timescale 
icoi ~ 400 (0.9gcm-Vcr) Myr > tacc,p. 

3.4. Velocity Instability and Ejection of Planets 

Once the cohort of Neptune-mass oligarchs consumes 
-^1/2 the mass of the parent disk, they scatter one another 
onto highly elliptical and inclined orbits (G04, see their 
equation [111]; Kenyon and Bromley, 2006). This veloc- 
ity instability occurs because damping of planetary random 
velocities by dynamical friction with the disk can no longer 
compete with excitation by neighboring, crowded oligarchs. 

The epoch of large planetary eccentricities lasts until 
enough oligarchs are ejected from the system. We can esti- 
mate the ejection time by following the same reasoning that 
led to (14). Replace vk with the system escape velocity 
Wosc,sys ~ ^a, and replace Sp with the surface density of 
oUgarchs ~Mp/a^ (see equation [13]). Then solve for 



The coefficient of 0.1 is attributed to more careful ac- 
counting of encounter geometries; equation (15) gives ejec- 
tion times similar to those found in numerical simulations 
(G04). Neptune-mass oligarchs at a « 20 AU kick their ex- 
cess brethren out over feject ^ 600 Myr. Removal is faster 
if excess oligarchs are passed inward to Jupiter and Saturn. 

Oligarchs moving on eccentric orbits likely traverse dis- 
tances beyond 30 AU and stir KBOs. We expect more mem- 
bers are added to the Scattered KBO disk during this stage. 

We have painted a picture of dynamically hot oligarchs 
similar to that drawn by Thommes et al. (1999; see also 
Tsiganis et al., 2005), who hypothesize that Neptune and 
Uranus form as oligarchs situated between the cores of 
Jupiter and Saturn at 5-10 AU. The nascent ice giants are 
scattered outward onto eccentric orbits once the gas giant 
cores amass their envelopes. While Neptune and Uranus 
reside on eccentric orbits, they can stir KBOs in much the 
same way as we have described above (Thommes et al., 
2002). Despite the similarity of implications for the stirring 
of KBOs, the underlying motivation of the cosmogony pro- 
posed by Thommes et al. (1999) is the belief that Neptune- 
mass bodies do not form readily at distances of ~30 AU. 
Recent work highlighting the importance of inelastic col- 
lisions amongst very small bodies challenges this belief 
(G04; see §3.1). 

3.5. Dynamical Friction Cooling of Surviving Planets 

Planetary oligarchs that survive ejection — i.e., Uranus 
and Neptune — have their e's and i's restored to small values 
by dynamical friction with the renmant disk (comprising 
predominantly small KBOs of surface density a and veloc- 
ity dispersion u) over time 

_ Vp pRp f Vp Y nf.^ 

Idf.cool — — ~ — P7 , (lo) 

Vp ail \Vcsc,pJ 

where Rp, Vesc,p, and Vp ^ u arc the planet's radius, sur- 
face escape velocity, and random velocity, respectively. For 
Vp = ria/2 (planetary eccentricity Cp ^ 0.5), a = 25 AU, 
Rp = 25000 km, Vesc,p = 24kms~^, and cr = Sp = 
0.9gcm~^ (since the velocity instability occurs when the 
surface density of oligarchs equals that of the parent disk; 
§3.4), we find tdf.cooi ~ 20 Myr. 

While Neptune's orbit is eccentric, the planet might re- 
peatedly invade the Kuiper belt at a w 40-45 AU and stir 
KBOs. Neptune would have its orbit circularized by trans- 
ferring energy to both small and large KBOs. Unlike small 
KBOs, large ones cannot shed this energy because they cool 
too inefficiently by inelastic collisions (see the end of §3.3). 
Insert (16) into (14) and set Sp = ct to estimate the ran- 
dom velocity to which large KBOs are excited by a cooling 
Neptune: 

t^K ~ Vp • (17) 

Thus large KBOs are stirred to the same random velocity 
that Neptune had when the latter began to cool, regardless 
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of the numerical value of idf.cooi- Large KBOs effectively 
record the eccentricity of Neptune just prior to its cooling 
phase. Final eccentricities ck might range from ~0.1 to 
nearly 1. During this phase, the population of the Scat- 
tered KBO disk would increase, perhaps dramatically so. 
If all large KBOs are stirred to ck ^ 0.1, new large KBOs 
must coagulate afterwards from the remnant disk of small, 
dynamically cold bodies to re-constitute the cold Classical 
disk. Cold Classicals might therefore post-date hot KBOs. 



disk of small bodies is unknown (§3.6.4). 

In addition to scouring the disk, Neptune's migration 
has been proposed to sculpt the disk in other ways — ^by 
capturing bodies into mean-motion resonances (§3.6.1), re- 
distributing the Classical disk by resonance capture and re- 
lease (§3.6.2), and deflecting objects onto Scattered orbits 
(§3.6.3). We critically examine these proposals below. 

3 .6. 1 . Capture and Excitation of Resonant KBOs 



3.6. Planetary Migration 

Having seen a few of its siblings evicted, and having set- 
tled onto a near-circular, flattened orbit, Neptune remains 
immersed in a disk of small bodies. The total mass of the 
disk is still a few times that of the planet because the prior 
velocity instability occurred when the surface density of oli- 
garchs was comparable to that of the disk. By continuing to 
scatter small bodies, Neptune migrates: Its semi-major axis 
changes while its eccentricity is kept small by dynamical 
friction. Absent other planets, migration would be Sunward 
on average. Planetesimals repeatedly scattered by Neptune 
would exchange angular momentum with the planet in a 
random-walk fashion. Upon gaining specific angular mo- 
mentum ~(\/2 — 1)^0^, where fl and a are appropriate 
to Neptune's orbit, a planetesimal initially near Neptune 
would finally escape. Having lost angular momentum to the 
ejected planetesimal, Neptune would migrate inward. (A 
single planet can still migrate outward if it scatters material 
having predominantly higher specific angular momentum. 
Gomes et al. (2004) achieve this situation by embedding 
Neptune in a disk whose mass is at least lOOM^ and is 
weighted towards large distances [aa^ oc a]; see also the 
chapter by Morbidelli et al.). 

Other planets complicate this process. Numerical simu- 
lations by Fernandez andip (1984) and Hahn andMalhotra 
(1999) incorporating all 4 giant planets reveal that planetes- 
imals that originate near Neptune are more likely ejected by 
Jupiter. Over the course of their random walks, planetesi- 
mals lose angular momentum to Neptune and thereby cross 
Jupiter's orbit. Jupiter summarily ejects them (see equation 
[15] and related discussion). Thus, on average, Neptune 
gains angular momentum and migrates outward, as do Sat- 
urn and Uranus, while Jupiter's orbit shrinks. 

An outward bound Neptune passes objects to the interior 
planets for eventual ejection and seeding of the Oort Cloud. 
We refer to this process as "scouring" the trans-Neptunian 
disk. Scouring and migration go hand in hand; the fraction 
by which Neptune's semi-major axis increases is of order 
the fraction that the disk mass is scoured. Scouring is likely 
a key part of the solution to the clean-up (a.k.a. missing- 
mass) problem. If clean-up is not achieved by the end of 
Neptune's migration, one must explain how to transport 
the bulk of the trans-Neptunian disk to other locales while 
keeping Neptune in place {Gomes et al., 2004). Scouring 
has only been treated in collisionless N-body simulations. 
How scouring and migration proceed in a highly coUisional 



As Neptune migrates outward, its exterior mean-motion 
resonances (MMRs) sweep across trans-Neptunian space. 
Provided the migration is sufficiently slow and smooth, 
MMRs may trap KBOs and amplify their orbital eccen- 
tricities and, to a lesser extent, their inclinations. The ec- 
centric orbits of Pluto and the Plutinos — objects which all 
inhabit Neptune's 3:2 resonance — may have resulted from 
resonance capture and excitation by a migrating Neptune 
(Malhotra, 1993, 1995; Jewitt and Luu, 2000). The ob- 
served occupation of other low-order resonances — e.g., the 
4:3, 5:3, and 2:1 MMRs— by KBOs on eccentric orbits (see 
Fig. 2 and Table 1) further support the migration hypothe- 
sis (C03). In this section, we review the basic mechanism 
of resonant excitation of eccentricity, examine how the mi- 
gration hypothesis must change in Ught of the unexpected 
occupation of high-order(e.g., the 7:4, 5:2, and 3:1) MMRs, 
and discuss how m: 1 resonances serve as speedometers for 
Neptune's migration. 

Consider the interaction between a test particle (KBO) 
and a planet on an expanding circular orbit. In a frame of 
reference centered on the Sun and rotating with the planet's 
angular velocity fip(t), the particle's Hamiltonian is 

n = £ -np{t)c-'R{t), (18) 

where £ = -GMQ/2a, C = [GMQa{l - e^)^!'^, and U 
is the disturbing potential due to the planet (these quanti- 
ties should be expressed in canonical coordinates). From 
Hamiltonian mechanics, dH/dt = dH/dt = —SlpC — 
dTZ/dt. Therefore 



(19) 



where e = {dTl/dt - dH/dt)/{d£/dt). We re-write (19): 



d^ 
dt 



(l-e^) 



1/2 



(1 - e")i/2 - n/np{i - e) 



da 



(20) 

where Q. is the particle's angular frequency. 

For a particle trapped in m:n resonance (where m and 
n are positive, relatively prime integers), a, e, and the res- 
onance angle change little over the particle's orbital period. 
If the synodic period is not much longer than the orbital 
period, we may average the Hamiltonian over the former 
(we may do this by choosing appropriate terms in the ex- 
pansion of TZ). This yields Q./Q.p{l — e) = n/m. For a 
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particle in resonance, |e| <C 1. By change of variable to 
a; = (1 - e^)^/^, equation (20) integrates to 



(1 - e2)i/2 _ 



a = constant , 



(21) 



which relates changes in a to changes in e for any 
resonance — exterior m > n, interior m < n, or Trojan 
m — n. In the case of a planet that migrates towards 
a particle in exterior resonance, a increases to maintain 
resonant lock {Goldreich, 1965; Peale, 1986). Then by 
(21), e also tends to increase, towards a maximum value 
[1 — (n/w)^]^/^. Particles inhabiting either an exterior or 
interior resonance have their eccentricities ampUfied from 
because they are perturbed by a force pattern whose angu- 
lar speed Qp does not equal their orbital angular speed f2. 
Particles receive energy and angular momentum from the 
planet in a ratio that cannot maintain circularity of orbits. 

Among observed 2:1 Resonant KBOs, max(e) ~ 0.38 
(Fig. 2). If 2:1 Resonant KBOs had their eccentricities am- 
plified purely by migration, they must have migrated by 
Aa « 13 AU (equation [21]). Neptune must have migrated 
correspondingly by Aop « 8 AU. This is an upper bound 
on Aflp because it does not account for non-zero initial ec- 
centricities prior to capture. 

In early simulations {Malhotm, 1993, 1995) of reso- 
nance capture by a migrating Neptune, resonances swept 
across KBOs having initially small e's and i's. These mod- 
els predicted that if Neptune's orbit expanded by Aop « 
8 AU, low-order resonances such as the 4:3, 3:2, 5:3, and 
2:1 MMRs would be occupied by objects having 0.1 < e < 
0.4andi < 10°. Eccentric KBOs indeed inhabit these reso- 
nances (Fig. 2). Two observations were not anticipated: (1) 
Resonant KBOs are inclined by up to i ~ 30°, and (2) high- 
order resonances — e.g., the 5:2, 7:4, and 3:1 — enjoy occu- 
pation. These observations suggest that Neptune's MMRs 
swept across not only initially dynamically cold objects, but 
also initially hot ones: The belt was pre-heated. For ex- 
ample, to capture KBOs into the 5:2 MMR, pre-heated ec- 
centricities must be > 0.1 (C03). Neptune-sized perturbers 
(§3.3-3.5) might have provided the requisite pre-heating in 
e and i. 

To understand why capture into high-order resonances 
favors particles having larger initial e, recognize that cap- 
ture is only possible if, over the time the planet takes 
to migrate across the maximum possible libration width 
max(^aiib), the particle completes at least 1 Ubration: 



max((5aiib) 



lib > 



(22) 



where Tub is the libration period {Dermott et al., 1988). 
Otherwise, the particle would hardly feel the resonant per- 
turbation as the planet races towards it. Since max(feiib) ~ 
('rorb/Tiib)ap and Tub - T,,,,{MQe-\'^-"\ / M,,^' 
{Murray and Dermott, 1999), where Torb is the orbital pe- 
riod of the particle and Mp is the mass of the planet, we 
re- write equation (22) as 



rpl 

-'lib 
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TorbTniig 

where the migration timescale Tmig = ap/|ap|. The higher 
the order |m — n| of the resonance, the greater e must be to 
satisfy (23) (C03; Hahn and Malhotra, 2005). 

Asymmetric (m:l) resonances afford a way to estimate 
the migration timescale observationally. An asymmetric 
MMR furnishes multiple islands of libration. At the fixed 
point of each island, a particle's direct acceleration by Nep- 
tune balances its indirect acceleration by the Sun due to the 
Sun's reflex motion (Pan and Sari, 2004; Murray-Clay and 
Chiang, 2005, hereafter MC05). The multiplicity of islands 
translates into a multiplicity of orbital longitudes, measured 
relative to Neptune's, where resonant KBOs cluster on the 
sky. The pattern of clustering varies systematically with 
migration speed at the time of capture (Chiang and Jordan, 
2002). For example, when migration is fast — occurring on 
timescales Tmig ^ 20 Myr — objects are caught into 2: 1 res- 
onance such that more appear at longitudes trailing, rather 
than leading, Neptune's. The degree of asymmetry can be 
as large as 300%. When migration is slow, the distribu- 
tion of captured 2:1 objects is symmetric about the Sun- 
Neptune line. The preference for trailing versus leading 
longitudes arises from migration-induced shifts in the sta- 
ble and unstable equiUbria of the resonant potential. Shifts 
in the equilibrium values of the resonance angle are given 
in radians by equation (23) and are analogous to the shift in 
the equilibrium position of a spring in a gravitational field 
(MC05). The observation that trailing 2: 1 KBOs do not out- 
number leading ones constrains Tmig > 20 Myr with nearly 
3ct confidence (MC05). This measurement accords with nu- 
merical simulations of the migration process itself by Hahn 
and Malhotra (1999) and by Gomes et al. (2004, see their 
fig. 10); in these simulations, Tmig ^ 40 Myr. 

3.6.2. Stochastic Migration and Resonance Retainment 

Finite sizes of planetesimals render planetary migra- 
tion stochastic ("noisy"). The numbers of high and low- 
momentum objects that Neptune encounters over fixed 
time intervals fluctuate randomly. These fluctuations spo- 
radically hasten and slow — and might occasionally even 
reverse — the planet's migration. Apportioning a fixed disk 
mass to larger (fewer) planetesimals generates more noise. 
Extreme noise defeats resonance capture. Therefore the 
existence of Resonant KBOs — which we take to imply 
capture efficiencies of order unity — sets an upper limit 
on the sizes of planetesimals (small bodies) comprising 
the bulk of the mass of the disk. Murray-Clay and Chi- 
ang (2006, hereafter MC06) estimate this upper limit to be 
Smax ~ C(IOO) km; a shortened derivation of their result 
reads as follows. 

For a given planetesimal size, most noise is generated 
per unit mass disk by planetesimals having sub-Hill (u < 
wh,p = ^Rp/ct) velocity dispersions and semi-major axes 
displaced ±i?H,p from the planet's (MC06). A single such 



Mp e 



< 1 



(23) 
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planetesimal of mass fi, after undergoing a close encounter 
with the planet, changes the planet's semi-major axis by 
Aoi ~ ±(/u/Mp)i?H,p. The planet encounters such plan- 
etesimals at a rate N ~ aRj^ pflp/iJ,. Over the duration 
of migration ~(Aap/ap)Tniig, the planet's semi-major axis 
random walks away from its nominal (zero-noise) value by 
Attrnd ~ ±(A^Aaprmig/ap)i/2|Aai|. The Ubration am- 
pUtude in a of any resonant KBO increases by about this 
same | Aamd | • Then stochasticity does not defeat resonance 
capture if |Aarnd| < max((5aiib); that is, if 

< (MA ( pR.ea, \ 
'-[m^J [a^lpT^,,AaJ " ' ^^^^ 

which evaluates to s < O(100)km for Up = 30 AU, 
Aap = 8AU, Tmig « 40Myr, a = 0.2 gcm'^, and 
e = 0.2. 

The above constraint on size applies to those planetesi- 
mals that comprise the bulk of the disk mass. Noise is also 
introduced by especially large objects that constitute a small 
fraction of the disk mass. The latter source of noise has been 
invoked to explain the curious near-coincidence between 
the edge of the Classical disk (a = 47 AU) and Neptune's 2: 1 
resonance (a = 47.8 AU). Levison and Morbidelli (2003) 
suggest that the sweeping 2:1 MMR captures KBOs only to 
release them en route because of close encounters between 
Neptune and objects having ^10 times the mass of Pluto 
("super-Plutos"). Dynamically cold KBOs, assumed to co- 
agulate wholly inside 35 AU (§3.1.2), are thereby combed 
outwards to fill the space interior to the final location of the 
2: 1 MMR. Why the super-Plutos that are invoked to gener- 
ate stochasticity have not been detected by wide-field sur- 
veys is unclear {Morbidelli et al, 2002). The scenario fur- 
ther requires that ^3M© be trapped within the 2: 1 MMR 
so that a secular resonance maintains a population of 2:1 
resonant KBOs on low-e orbits during transport. 

3.6.3. Contribution of Migration to Scattered KBOs 

Neptune migrates by scattering planetesimals. What 
fraction of these still reside today in the Scattered belt? Do 
hot Classicals (having i > 5°) owe their excitation to a mi- 
gratory Neptune? Many Scattered and hot Classical KBOs 
observed today have q > 37 AU. This fact is difficult to ex- 
plain by appealing to perturbers that reside entirely inside 
30 AU. Insofar as a close encounter between a perturber 
and a particle can be modelled as a discontinous change in 
the particle's velocity at fixed position, the particle (assum- 
ing it remains bound to the Sun) tends to return to the same 
location at which it underwent the encounter. 

Gomes (2003ab) proposes that despite this difficulty, ob- 
jects scattered by Neptune during its migration from ^20 
to 30 AU can evolve into today's Scattered and hot Clas- 
sical KBOs by having their perihelia raised by a variety of 
sweeping secular resonances (SRs; see §2.1). As the outer 
planets migrate, SRs sweep across trans-Neptunian space. 



After having its e and i amplified by close encounters with 
Neptune, a planetesimal may be swept over by an SR. Un- 
like MMRs, SRs cannot alter particle semi-major axes and 
therefore do not permanently trap particles. However, a par- 
ticle that is swept over by an apsidal-type SR can have its 
eccentricity increased or decreased. A particle swept over 
by a secular resonance is analogous to an ideal spring of 
natural frequency uq, driven by a force whose time- variable 
frequency u){t) sweeps past uq. Sweeping uj past wq can 
increase or decrease the amphtude of the spring's free os- 
cillation (the component of the spring's displacement that 
varies with frequency loq), depending on the relative phas- 
ing between driver and spring near the moment of resonance 
crossing when uj ojq. 

Lowering e at fixed a raises q. Gomes (2003ab) and 
Gomes et al. (2005) find in numerical simulations of plan- 
etary migration that Neptune-scattered planetesimals origi- 
nating on orbits inside 28 AU can have their periheha raised 
up to 69 AU by a combination of sweeping SRs, MMRs, 
and Kozai-type resonances (which are a kind of SR). In ad- 
dition to offering an explanation for the origin of high-g, 
high-z KBOs, this scenario also suggests a framework for 
understanding differences in physical properties between 
dynamical classes. Compared to Classical KBOs, which 
are held to coagulate and evolve largely in situ. Scattered 
KBOs originate from smaller heliocentric distances d. To 
the (unquantified) extents that coagulation rates and chemi- 
cal environments vary from ~ 20 to 50 AU, we can hope 
to understand why a large dispersion in i — which in the pro- 
posed scenario reflects a large dispersion in birth distance 
d — implies a large dispersion in color/size. 

The main difficulty with this perihelion-raising mecha- 
nism is its low efficiency: Only ~0.1% of all objects that 
undergo close encounters with a migratory Neptune have 
their perihelia raised to avoid further close encounters over 
the age of the solar system (Gomes, 2003ab). Based on 
this mechanism alone, a disk weighing ~5OM0 prior to 
migration would have ^O.G5M© deposited into the Scat- 
tered and hot Classical belts for long-term storage. But 
only ~l-2% of this mass would be in bodies having sizes 
R > 100 km (/Tenyon and Luu, 1998, 1999; Kenyan, 2002; 
§3.1). Therefore this scenario predicts that Scattered and 
hot Classical KBOs having R > 100 km would weigh, in 
total, ^1Q~^M^ — about 50-150 times below what is ob- 
served (§2.2). This discrepancy is missed by analyses which 
neglect consideration of the KBO size distribution. A sec- 
ondary concern is that current numerical simulations of this 
mechanism account for the gravitational effects of disk par- 
ticles on planets but not on other disk particles. Proper 
calculation of the locations of secular resonances requires, 
however, a full accounting of the mass distribution. 

Given the low efficiency of the mechanism, we submit 
that the high-(j orbits of hot Classical and Scattered KBOs 
did not arise from Neptune's migration. Instead, these or- 
bits may have been generated by Neptune-mass oUgarchs 
whose trajectories passed through the Kuiper belt. While 
a numerical simulation is necessary to test this hypothesis. 
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our order-of-magnitude estimates (§§3.3-3.5) for the degree 
to which oligarchs stir the belt by simple close encounters 
are encouraging. No simulation has yet been performed in 
which the Kuiper belt is directly perturbed by a mass as 
large as Neptune's for a time as long as tdf ,cooi ~ 20 Myr. 
Differences in physical properties between Classical and 
Scattered/Resonant KBOs might still be explained along the 
same Unes as described above: Scattered/Resonant KBOs 
were displaced by large distances from their coagulation 
zones and so might be expected to exhibit a large disper- 
sion in color and size, while Classical KBOs were not so 
displaced. Even if all KBOs having R > 100 km were 
heated to large e or i by planetary oligarchs, the cold Clas- 
sical disk might have re-generated itself in a second wave 
of coagulation from a coUisional disk of small bodies. 

3.6.4. Problems Regarding Migration 

The analyses of migration cited above share a common 
shortcoming: They assume that planetesimals are collision- 
less. But coagulation studies (§3.1) indicate that much of 
the primordial mass remains locked in small bodies for 
which collision times threaten to be shorter than the du- 
ration of planetary migration. By (3), planetesimals hav- 
ing sizes <C 1 km in a minimum-mass disk have colUsion 
times <§; 20 Myr. How Neptune's migration unfolds when 
most of the disk comprises highly coUisional bodies has not 
been well explored. Neptune may open a gap in the disk 
(in the same way that moons open gaps in coUisional plan- 
etary rings) and the planet's migration may be tied to how 
the disk spreads by colUsional diffusion {Goldreich et al, 
2004b). 

How does the Classical belt shed 99% of its primordial 
mass? Situated at 40^7 AU, it may be too distant for Nep- 
tune to scour directly. Perhaps the small bodies of the Clas- 
sical belt are first transported inwards, either by gas drag or 
coUisional diffusion, and subsequently scoured. Clean-up 
and migration are intertwined, but the processes are often 
not discussed together (but see Gomes et al, 2004). 

Are there alternatives to migration for the capture of 
Resonant KBOs? Perhaps Resonant KBOs are captured 
as Neptune's orbit cools by dynamical friction (§3.5). Be- 
fore capture, many belt members would already be stirred to 
large e and i, not only by unstable oligarchs (§3.4), but also 
by Neptune while it cools. Cooling accelerates as it pro- 
ceeds (equation [16]). A rapid change in the planet's semi- 
major axis towards the end of cooling might trap KBOs 
into resonance by serendipity. Just after Neptune's semi- 
major axis changes, objects having orbital elements (includ- 
ing longitudes) suitable for libration would be trapped. This 
speculative "freeze-in" mechanism might be too inefficient, 
since it requires that the fraction of phase-space volume oc- 
cupied by resonances equal the fraction of KBOs that are 
Resonant. Taken at face value, observations suggest the lat- 
ter fraction is not much smaller than order unity (§2.1). 

3.7. Stellar Encounters 



A passing star may have emplaced Sedna onto its high- 
perihelion orbit. For the last t ~ 4 Gyr, solar-mass stars 
in the Solar neighborhood have had an average density 
n* ~ 0.04 stars pc~^ and a velocity dispersion (wj)^/^ ~ 
30kms~^. If we assume that the Sun once resided within 
a "typical" open cluster, then n* ~ 4 stars pc~^ and 
{vlY'"^ ~ Ikms"^ over t ~ 200 Myr. Over t > 
200 Myr, open clusters dissolve by encounters with molec- 
ular clouds (Binney and Tremaine, 1987). The number 
of stars that fly by the Sun within a distance 5* large 
enough that gravitational focussing is negligible (g, > 
GMq/{vI) ~ 900AUfor (?;2)i/2 ^ 1 km g-^) increases 
as J* ni,{v'^y^^dt. Therefore fly-bys during the current 
low-density era outnumber those during the cluster era by 
a factor of ~6. Nonetheless, intra-cluster encounters can 
be more effective at perturbing KBO trajectories because 
encounter velocities are 30 x lower. 

Fernandez and Brunini (2000) simulate the formation 
of the Oort cloud within an open cluster having param- 
eters similar to those cited above. They find that pass- 
ing stars create an "inner Oort cloud" of objects having 
35 < (7(AU) < 1000, 300 < a(AU) < 10^ (e) ~ 0.8, and 
^^2^1/2 ^ 2. Sedna may be the first discovered member of 
this inner Oort cloud {Brown et al, 2004). Such objects co- 
agulate in the vicinity of the giant planets and are scattered 
first by them. Since a scattering event changes velocities 
more effectively than it does positions, objects' perihelia 
remain at heliocentric distances of ~5-30 AU while aphe- 
lia diffuse outward. Aphelia grow so distant that objects 
are scattered next by cluster stars. These stars raise objects' 
periheUa beyond the reach of the giant planets. 

We confirm the ability of cluster stars to raise the perihe- 
lion of Sedna with an order-of-magnitude calculation. Dur- 
ing the open cluster phase, the number of stars that pass 
within distance of the Sun is 



V4OOOAU/ 



4pc-3 lkms-1 200 Myr 



(25) 

A star of mass A/, having perihelion distance (7* much 
greater than a planetesimal's aphelion distance (Q ~ 2a) 
perturbs that object's specific angular momentum by 



5h = ±C- 



(26) 



where the numerical coefficient C depends on the encounter 
geometry (Yabushita, 1972). We can derive the form of (26) 
by noting that ~ Q6v, where 5v is the perturbation to the 
object's velocity relative to the Sun. We write Sv as the tidal 
acceleration GM^Q/q^ induced by the star, multiplied by 
the duration of the encounter, to arrive at (26). 

For highly eccentric orbits 6q = hSh/ (GMq), whence 



5q 
Q 



±C 



Mr. 



/2GM0 



1/2 



(27) 
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For = Mq, = 4000 AU, C w 6 (see equa- 
tion [3.17] of Yabushita [1972]), {vlf/"^ = lkms-\ and 
pre-encounter values of q = 35 AU and a = 600 AU, 
5q/q ~ ±1. Thus, Sedna's perihelion could have dou- 
bled to near its current value, q « 76 AU, by a single 
slow-moving cluster star. Multiple encounters at larger 
cause q to random walk and change its value less effec- 
tively: {{5q)Y'^ oc {N^f'^z'' oc q-\ 

Had we performed this calculation for parameters ap- 
propriate to the present-day stellar environment, we would 
have found 5q/q « ±0.2. The reduction in efficacy is due 
to the larger {v1) today. 

The cluster properties cited above are averaged over a 
half-light radius of 2 pc (Binney and Tremaine, 1987). For 
comparison, the Hyades cluster has 4x lower n,, 3x lower 
(wj)^/^, and 6x longer lifetime t {Binney and Merrifield, 
1998; Perryman et ai, 1998); the Hyades therefore gener- 
ates 2x fewer encounters than does our canonical cluster. 
Younger clusters like the Orion Trapezium maintain 15 x 
higher n* and similar (wj)^/^ over 200x shorter t (Hil- 
lenbrand and Hartmann, 1998), and therefore yield even 
fewer encounters. Scenarios that invoke stellar encounters 
for which <C 1000 AU to explain such features as the 
edge of the Classical belt require that the Sun have resided 
in a cluster having atypical properties, i.e., dissimilar from 
those of the Orion Trapezium, the Hyades, and all open 
clusters documcntedby Binney and Merrifield (1998). That 
parent bodies in extra-solar debris disks also do not extend 
beyond ^40-100 AU (§3.1.2) argues against explanations 
that rely on unusually dense environments. 



In the theory of oligarchic planet formation (e.g., G04), 
each annulus is of order 5i?H in radial width; the number of 
QR-sized oligarchs that can be fitted into the tadpole Ubra- 
tion region is 



Trojan ' 



(8Mn/3Mq) 



1/2 



ON 



5Rk 



20, 



(28) 



attractively close to the number of QR-sized Neptune Tro- 
jans inferred to exist today (§2.2.3). The numerator in (28) 
equals the maximum width of the 1 : 1 MMR, an ~ 30 AU 
is Neptune's current semi-major axis, i?H = i?/a is the 
Trojan's Hill radius, and i? « 90 km is the radius of QR. 

The input parameters of the coagulation model are the 
surface density cr and sizes s of small bodies in 1:1 reso- 
nance. Big bodies grow by consuming small bodies, but 
growth is limited because small bodies diffuse out of reso- 
nance by colliding with other small bodies. The time for a 
small body to random walk out of the Trojan sub-disk is 



ps_ 

an 



(29) 



The term in square brackets follows from noting that a 
small body shifts its orbital guiding center by of order its 
epicyclic amplitude ~±u/r2 every time it collides with an- 
other small body in an optically thin disk. To escape res- 
onance, the small body must random walk the maximum 
libration width. We equate tesc to the growth time of a big 
body ^acc (equation [4]) to solve for the maximum size to 
which a large body coagulates: 



3.8. Coagulation of Neptune IVojans 

Planetesimal collisions that occur near Neptune's La- 
grange points insert debris into 1 : 1 resonance. This debris 
can coagulate into larger bodies. The problem of accretion 
in the Trojan resonance is akin to the standard problem of 
planet formation, transplanted from a star-centered disk to a 
disk centered on the Lagrange point. As with other kinds of 
transplant operations, there are compUcations: Additional 
timescales not present in the standard problem, such as the 
libration period Tub about the Lagrange point, require jug- 
gling. Chiang andLithwick (2005, hereafter CL05) account 
for these complications to conclude that QR-sized Trojans 
may form as miniature ohgarchs, each dominating its own 
tadpole- shaped annulus in the ancient Trojan sub-disk. Al- 
ternative formation scenarios for Trojans such as pull-down 
capture and direct collisional emplacement of QR-sized ob- 
jects into resonance are considered by CL05 and deemed 
unlikely. Also, the mechanism proposed by Morbidelli et 
al. (2005) to capture Jupiter Trojans cannot be applied to 
Neptune Trojans since Uranus and Neptune today he inside 
their 1:2 MMR and therefore could not have divergently 
migrated across it (Morbidelli, personal communication). 
We focus on in situ accretion, but acknowledge that a coUi- 
sionless capture scenario might still be feasible and even fa- 
vored by late-breaking data; see the end of this sub-section. 



R = i?final ~ 100 



\u/vuj V 



20 cm 



1/3 



km . (30) 



Our normalization of u/wh 



and (7 ~ 4 X 10 ^gcm ^ 



2 is derived from s ^ 20 cm 
^10 times the surface den- 
sity inferred in QR-sized objects today; we derive u/vuhy 
balancing gravitational stirring by big bodies with damping 
by inelastic colhsions between small bodies (CL05). For 
these parameter values, t 

esc tacc ~ 1 X 10^ yr. Unlike 
Neptune-sized oligarchs that may have been ejected out of 
the solar system (§3.4), all ~ 10-30 Trojan oligarchs in a 
single cloud should be present and eventually accounted for. 

As speculated by CL05, orbital inclinations of Trojans 
with respect to Neptune's orbit plane might be small; per- 
haps (i^)^/^ < 10°. A thin disk of Neptune Trojans would 
contrast with the thick disks occupied by Jupiter Trojans, 
main belt asteroids, and non-Classical KBOs, and would 
reflect a collisional, dissipative birth environment. Three 
other Neptune Trojans have since been announced after the 
discovery of QR, having inclinations of 1.4°, 25.1°, and 
5.3° (Sheppard and Trujillo, 2006). If a large fraction of 
Neptune Trojans have high i, we might look to the sec- 
ular resonance, unmodelled by CL05, to amplify inchna- 
tions. See also Tsiganis et al. (2005) who find that Neptune 
Trojans can be captured collisionlessly; the capture process 
is related to "freeze-in" as described in §3.6.4. 
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3.9. Collisional Comminution 

Over the last few billion years, sufficiently small and nu- 
merous bodies in the Kuiper belt suffer collisional attrition. 
As interpreted by Pan and Sari (2005, hereafter PS05), the 
break in the size distribution of KBOs at i? « 50 km as 
measured by Bernstein et al. (2004; §2.2.1) divides the col- 
lisional spectrum at small R from the primordial coagula- 
tion spectrum at large R. For the remainder of this subsec- 
tion, we do not distinguish between the various dynamical 
classes but instead analyze all KBOs together as a single 
group. At i? > -Rbreak. the size spectrum dN/dR oc R~^°, 
where dN is the number of objects per unit face-on area 
of the belt having sizes between R and R + dR (the dif- 
ferential surface number density). The slope (?o ~ 5 (see 
§2.2.1 for more precise values) presumably represents the 
unadulterated outcome of coagulation. Bodies at this large- 
R end of the spectrum are insufficiently numerous to colhde 
amongst themselves and undergo attrition. At i? < i?break, 
dN/dR oc i?^', where q derives from a quasi-steady col- 
lisional cascade (Dohnanyi, 1969; PSQ5). By definition of 
-Rbreak, the time for a body of radius i?break to be catas- 
trophically dispersed equals the time elapsed: 



We assume (and can check afterwards) that i?. 



t, 



(31) 



where T^R^^^i-ik the collision cross-section and Nproj is 
the surface number density of projectiles that are just large 
enough to disperse -Rbreak-sized targets (catastrophic dis- 
persal implies that the mass of the largest post-collision 
fragment is no greater than half the mass of the original 
target and that collision fragments disperse without gravi- 
tational reassembly). This expression is valid for the same 
assumptions underlying equation (3), i.e., for today's dy- 
namically hot belt. 

We proceed to estimate i?broak given the parameters 
of the present-day Kuiper belt. For R > i?brcak, N = 
No{R/ Rq)^^''" , where N is the surface number density 
of objects having sizes between R and 2R. We estimate 
that for fiducial radius Rq = 100 km, iVo ~ 20 AU"^ at 
a i=a 43 AU. The minimum radius .Rproj of the projectile 
that can catastrophically disperse a target of radius inbreak 
is given by 



- ff3 iP 

2 proj ^rel 



-^breakQ* 



where 



Q* 



(32) 



(33) 



is the collisional specific energy (Greenberg et al., 1978; 
Fujiwara et al., 1989) and t^rei is the relative collision ve- 
locity. Since for R < .Rbreak as much mass is ground into 
every logarithmic interval in as is ground out (e.g., PSQ5), 
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Combining the above relations yields 
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(36) 



where 2;i = (6 -I- y)/[^y + (6 + y){qo - 3)] and Z2 = 
5/[5?y+(6 + j/)((Jo — 3)]. For targets held together by self- 
gravity, Q* « 3w^g(,/10 and y = 2. If we insert these 
values into (36), together with v^ei = Ikrns"^, go = 5, 
n = 27r/(300yr), and t = 3 x 10^ yr, we find that 
-Rbreak ~ 0.4i?o « 40 km, in good agreement with the ob- 
served break in the luminosity function (Fig. 5; PS05). The 
small- i? end of the KBO size spectrum as observed today 
reflects the catastrophic comminution of bodies that derive 
their strength from self-gravity ("rubble piles"). Further- 
more, the Kuiper belt has been dynamicaUy hot for the last 
few bilhon years (PS05). 

4. DIRECTIONS FOR FUTURE WORK 

1. Collisional vs. Collisionless: Most explorations of 
planetary migration and of how the Kuiper belt was stirred 
utilize colhsionless gravitational simulations. But the over- 
whelming bulk of the primordial mass may have resided 
in small, collisional bodies. Simultaneously accounting for 
collisions and gravity might revolutionize our understand- 
ing of the clean-up (a.k.a. missing-mass) problem. Insights 
from the study of planetary rings will be helpful. 

2. Classical KBO Colors vs. Heliocentric Distance: Do 
Classical KBOs exhibit a trend in color from neutral to red 
with increasing heliocentric distance dl The two neutral 
Classicals at d ~ 38 AU, contrasted with the predomi- 
nantly red Classicals at d « 42 AU, suggest the answer is 
yes (Figs. 3 and 4). Confirmation would support ideas that 
Classicals coagulated in situ, and that neutrally colored Res- 
onant/Scattered KBOs coagulated from small d and were 
transported outwards. We must also ask why trends in color 
with birth distance d would exist in the first place. 

3. Formation of the Scattered Belt by Neptune-Mass Oli- 
garchs: We argue that Neptune's migration and the con- 
comitant sweeping of secular resonances do not populate 
the Scattered and hot Classical belts with enough objects to 
explain observations. When account is made of the primor- 
dial size distribution of planetesimals — a distribution that 
should be preserved today at large sizes (§1; §3.9) — ^the ex- 
pected population of Scattered / hot Classical objects hav- 
ing sizes above 100 km is less than that observed by a fac- 
tor of 50-150. We propose instead that planetesimals were 
deflected onto Scattered / hot Classical orbits by simple 
close encounters with marauding Neptune-mass oligarchs 
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that have since been ejected from the solar system, and by 
Neptune while its orbit circularized by dynamical friction. 
These contentions are supported by order-of-magnitude es- 
timates but require numerical simulations to verify. 

4. Kuiper Cliff: Why do planetesimal disks have sharp 
outer edges? 

5. Binaries: Kuiper belt binaries might prove the most 
informative witnesses we have to the history of trans- 
Neptunian space. They hearken back to a primordially 
dense and cold disk in which colUsions and multiple-body 
encounters were orders of magnitude more frequent than 
they are today. Binary orbit properties must also reflect 
how the Kuiper belt was stirred as a whole. How binary 
inclinations, eccentricities, and component mass ratios are 
distributed, and how/why the incidence of binarity corre- 
lates with dynamical class are open issues for observer and 
theorist alike. 
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Fig. 1. — Orbital elements, time-averaged over 10 Myr, of 529 securely classified outer solar system objects as of 8 Oct 2005. Symbols 

represent dynamical classes: Centaurs (x), Resonant KBOs (O), Classical KBOs (•), and Scattered KBOs (o). Dashed vertical lines 
indicate occupied mean-motion resonances; in order of increasing heliocentric distance, these include the 1:1, 5:4, 4:3, 3:2, 5:3, 7:4, 9:5, 
2:1, 7:3, 5:2, and 3:1 (see Table 1). Solid curves trace loci of constant perihelion q = a(l — e). Especially large (2003 UB313, Pluto, 
2003 ELei, 2005 FYg; Brown et al, 2005ab) and dynamically unusual KBOs are labelled (2001 QR322 [Trojan; Chiang et al., 2003; 
Chiang and Lithwick, 2005], 2000 CRios [high q; Millis et al, 2002; Gladman et al, 2002], Sedna [high q; Brown et al, 2004]). For a 
zoomed-in view, see Figure 2. 
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